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ABSTRACT 
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I. Introduction 

Temperature induced symmetry-changing phase-transitions in quantum field theory 
are important ingredients in cosmological scenarios. The existence of the high-temperature 
phase transitions was suggested by Kirzhnits and Linde[l], and was shown quantitatively 
by Dolan and Jackiw[2], Weinberg [3], and Kirzhnits and Linde[4], by calculating the finite 
temperature one-loop effective potential. However, cosmological scenarios often rely on the 
detailed nature of the phase transition, namely whether it is of first or second order. A more 
precise determination of the critical temperature and the nature of the phase transition 
requires an analysis of higher-loop contributions even when the coupling constants in the 
theory are very small. Weinberg has argued, by using power counting, that the leading 
contributions at very high temperatures come from all those loops with superficial degree 
of divergence D> 1. This implies that in order to obtain more accurate information one 
must study infinite series of certain classes of multi-loop diagrams in perturbation theory. 
For example, in scalar theory the leading contributions come from the multi-loop 
graphs shown in Fig.l, which are called daisy and super-daisy. For this reason, Dolan 
and Jackiw, in their early paper, studied the effect of these graphs on the temperature 
dependent effective mass. 

Recently there has been much interest in the nature of the electroweak phase transi- 
tion due to the idea that the baryon asymmetry may be generated at the electroweak scale 
if the transition is of first order[5]. Several authors have calculated the high temperature 
effective potential in the standard model and in the theory, taking into account leading 
(and subleading) contributions from multi-loop diagrams, in order to obtain a correct form 
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of the high temperature effective potential [6, 13]. Some of the authors have calculated an 
"improved" one-loop effective potential in which the tree-level propagators are replaced by 
temperature dependent effective propagators, which were obtained by summing the dom- 
inant high temperature contributions from infinite-series of certain classes of self-energy 
graphs in perturbation theory. When one considers only the leading corrections to the 
effective propagators all results are in agreement with each other. However, there have 
been various disagreements when the subleading corrections to the effective propagators 
are included. On the other hand, the subleading contributions could be important in 
determining the nature of the phase transition. 

We find that in the improved one-loop calculations the difficulties arise due to the fact 
that the naive substitution of improved propagators in the one-loop effective potential is an 
ad-hoc approximation. One needs a self-consistent loop expansion of the effective potential 
in terms of the full propagator. Such a technique has been developed some time ago by 
Cornwall, Jackiw and Tomboulis (CJT) in their effective action formalism for composite 
operators [14]. One considers a generalization r(0, G) of the usual effective action r(0) 
which depends not only on (j){x) -a possible expectation value of the quantum field $(x)- 
but also on G{x,y) -a possible expectation value of the time-ordered product T^{x)^{y). 
The physical solutions satisfy stationary requirements 



STicf), G) 



(1.1) 




SG{x,y) 



= 0. 
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The conventional effective action r{(f)) is given by r{(f),G) at the solution Go{(f)) of 



In this formalism, it is possible to sum a large class of ordinary perturbation-series di- 
agrams that contribute to the effective action r(0), and the gap equation which determines 
the form of the propagator is obtained by a variational technique, as in (1.2). 

For translationally invariant solutions, we set = constant, take G{x, y) to be a 
function of only {x — y), and obtain the effective potential for the full propagator: 



The purpose of this paper is to understand the structure of the leading corrections to 
the perturbative one-loop finite temperature effective potential for the theory using 
CJT formalism in imaginary (Euclidean) time. We obtain a finite temperature effective 
potential, which is exact up to the order that includes all contributions from daisy and 
super-daisy graphs. Instead of dropping various finite and divergent terms, as has been 
done often in the recent literature, we carry out renormalization and then perform a high 
temperature expansion. We show explicitly that the effective potential must be calculated 
up to two-loop in order to generate all daisy and super-daisy graphs which appear in per- 
turbation theory[15]. Moreover, we find subtle cancellations of leading corrections between 
the improved one and two loop contributions. If it is indeed possible to determine the order 
of the electroweak phase transition by calculating the improved high temperature effective 
potential, our result implies that the improved two-loop contribution could play a crucial 
role. We plan to present our calculations in gauge theories in a future publication. 



(1.2) 



r(</>) = r(</>,Go(0)). 



(1.3) 




(1.4) 
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II. CJT Composite Operator Effective Potential 

In this section we shall review briefly the CJT formalism following Ref . [14] . To describe 
field theory at finite temperature T we shall use Euclidean time r, which is restricted to the 
interval 0<t</?=1/T. The Feynman rules are the same as those at zero temperature, 
except that the momentum space integral over the time component /C4 is replaced by a 
sum over discrete frequencies A;4 = 2TTn/j3 : 



where n is even (odd) for bosons (fermions). 

We introduce a partition function Zp{J^K) in the presence of the sources J and K 
defined by 



[We have set c = ^ = 1.] The ^-integration is functional and $ satisfies the periodic 
boundary conditions $(/3/2,x) = $(— /3/2,x). /($) is the classical Euclidean action, 
which may be written as 




(2.1) 




(2.2a) 





where Dq {x — y) is the free propagator 



Do\x -y) = -(□ + rr?)5\x - y) , 



(2.2c) 



and the interaction Lagrangian Li^t is at least cubic in $. 
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The effective action r/j(0, G) is obtained by a double Legendre transformation of 

W/3(0, G) = lnZp{(f)^ G) which is analogous to the free energy. We define 

5Wp{J,K) 



5J{x) 
5Wf3{J,K) 1 



5K{x,y) 2 
and eliminate J and K in favor of and G 



[G{x,y) + (j){x)(j){y)] , 



(2.3) 



^(0,G) = W^iJ,K) - J d'^x cl>{x)J{x) 



(2.4) 



It follows that 

^r((^,G) 

5(i){x) 



= J{x)- j d^yK{x,y)<j>{y) , (2.5) 

5G{x,y) 2 ^ ^ ' 

Since the physical processes correspond to vanishing sources J and Eqns.(2.5) and 
(2.6) provide the stationary requirement of (1. !)-(!. 2). r^(0, G) obtained in this way is 
the generating functional in </> for two-particle irreducible Green's functions expressed in 
terms of the full propagator G. 

In order to obtain a series expansion of r/3(0, G) we introduce the functional operator 
D-\cl)-x,y): 



The required series obtained by CJT is then 

r^(0, G) = Ici{4>) + ^TrLnDoG-^ + ^Tr[D-'G - 1] + rf \4>, G) , (2.8) 
where the (^-independent terms are chosen so that the overall normalization is consistent 

/3 



(2) 

with the conventional effective action r/3((/>) according to (1.3). The quantity r\'{(f)) is 
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computed as follows. In the classical action /($), shift the field $ by (p. Then /($ + (f)) 
contains terms cubic and higher in $ that define Iint{(i>'-i^) where the vertices depend on 
(f). (0) is given by all the two-particle irreducible (2PI) vacuum graphs in the theory 
with vertices given by Iint{4>; ^) and propagators set equal to G{x, y). 

From the stationary requirement (1.2) we obtain the gap equation for G: 

G-\.,y) = D-\.,y)-2^^. (2.9) 

When one is interested in translationally invariant solutions, the generalized effective 
potential V/3(0, G) can be obtained using (1.4) and (2.8). 

III. Theory 

A. Effective Potential ^^(0,(5) 
In this section we calculate the finite temperature effective potential for the single 
scalar field with A^"^ interaction. The Euclidean Lagrangian density is given by 

L=\{d,m''^) + \m''^'' + ^^' ■ (3.1) 
The propagator defined in (2.7) is 

D-^ {4>; X, y) = -(□ + + ^4^^)5\x - y) , (3.2) 

and the vertices of the shifted theory are given by 

i^int(<^;$) = ^0$' + ^$'. (3.3) 

In Fig. 2 the diagrams contributing to T^p\(j)^G) are shown up to three loops: each line 
represents the propagator G{x, y) and there are two kind of vertices. 
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In order to determine how to truncate the loop expansion so that aU daisy and super- 
daisy graphs of perturbation theory are included, we first observe that upon dropping F^^^ 
altogether, the gap equation (2.9) gives 



G-\x,y) = D-\x,y) , 



(3.4) 



and r{(f),D) is simply the ordinary one-loop effective potential. Therefore, a non-trivial G 
can be obtained only if we retain some of the graphs in \ Next we observe that among 
the graphs in Fig. 2 only the two-loop graph of 0(A) will include contributions from daisy 
and superdaisy graphs of ordinary perturbation theory. Therefore, we shall truncate the 
series at 0(A). This is known as Hartree-Fock approximation. The effective action is then 



By stationarizing F^ with respect to G we obtain the gap equation in the Hartree-Fock 
approximation: 



It is straightforward to show by iteration that Eqn.(3.6) generates all daisy and super-daisy 
graphs that contribute to the full propagator in ordinary perturbation theory. [Eqns.(3.5) 
and (3.6) have the same structure as those in the leading large N approximation since 
in both cases n-point functions are expressed in terms of one and two-point functions. 
However, as discussed in Ref.[14], in the large N approximation those daisy and super- 
daisy diagrams which are of 0{1/N) are dropped and therefore the coefficients of A in 
(3.5) and (3.6) become smaller by a factor three.] 





(3.6) 
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In order to obtain the effective potential V^(0, G) for translation-invariant field con- 
figurations, we define the Fourier-transformed propagators 

= ^J^^^. (3.7) 

where we have taken an Ansatz for G{k) using an effective mass M^. Since the correction 
to the gap equation in this approximation is given by XG{x,x)/2, can be taken to be 
momentum independent. The effective potential is then 

74) 



y (0, M) J-m^d^^ + ^ I -|A in[k^ + 



^[M^ - - ^(p^]G{x, x) + -G{x, x)G{x, x) . 
2 2 8 



(3.9) 



From (3.9) the stationary requirements are 

^ii^ = + + ^G{x, X)] = , (3.10) 

The conventional effective potential is obtained by evaluating Vg(0, M) at the solution 

M{(f)) of (3.11). It is composed of three terms: the classical (V^), one-loop (V^), and two- 
loop (V^^) contributions. 

V/3{4>, M{4>)) = + + V^^ , (3.12a) 

V' = \m^ct>^ + , (3.126) 

^'^/ (04 M/^' + M2(0)] , (3.12c) 

V^^ = -^G{x,x)G{x,x) , (3.12d) 
8 
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where 

M^{(P)^m^ + -(l)'^ + -G{x,x) , (3.12e) 

and now G{x, x) is given by 

B. Renormalizing the Effective Potential 

The expression of Va(0, M(0)) in (3.12) contains divergent integrals. Moreover, due 
to the fact that our approximation is non-perturbativcly self-consistent, reflecting the non- 
linearity of the full theory, M(0), the argument of Vg, is not well-defined due to the 
infinities in G{x,x). Therefore, we shall first obtain a well-defined, finite expression for 
M(0) by a renormalization. [In the rest of this paper M refers to the solution of (3.12e).] 

We define renormalized parameters mn and Xr as 

, m D 1 ^ , , 



l^ = i + i/2M, (3.13(.) 



where /i^2 are divergent integrals 



J (2^ ^2^ 



'^=y (2^^=ii^oo8^' ^^-^^^^ 

1 1 

, _ ] = lim :r In — ^ . (3.13(i) 

|k| A^ool67r2 p? ^ > 

fi is the renormalization scale and A is the ultraviolet momentum cut-off. [The same 

renormalization prescription has been used also in the large N approximation[17].] 

When the sum on n in A;4 is carried out as in Ref.[2], G{x, x) becomes 

^/c 1 



'^^'''^^-y (27r)3 2uu^ (3] (27r) 



r)3 Uk{exp[f3uk] - 1) (3.14) 
^G{M{ct^)) + h-M\ct))l2{ii) , 
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where tOk = [|kp + M^((/))]^/^ and G(M((/))) is the finite part of G{x,x) given by 

G(Mm . ^ m + - y ^^^^^^^^^j _ . (3.15) 

In the hmit T = 0, the first term in G{M{(j))) survives, but the second term vanishes. 

It is straightforward to see that (0) is finite and cut-off independent in terms of 
mR and \r 

M\cj>) = ±ml + + ^G(M(0)) . (3.16) 

[In this paper we shall choose the negative sign which allows spontaneous symmetry break- 
ing.] It is convenient for our later discussions to define the tree-level effective mass rn^{(f)) 

m\cj>) = -ml+^ct>' . (3.17) 

With this finite M^, we are ready to discuss the divergences in Va((/!), M). First, 
carrying out the sum on n in , we obtain the familiar one-loop finite temperature 
formula of Ref.[2], where the tree level effective mass m(0) is replaced by M(0) 
1 r d^k If d^k 



.ri^..^ i f d-'k If d'k , ,^ 



(3-18) 

M2 1 If d^k , r« n ^ w X 

= 64^[^" ^ - 2^ + ;5 y (2^ - - -T''^^^ ■ 

At T = the first term of survives and provides the zero-temperature one-loop contri- 
bution, and the second term vanishes. The last term is the divergence in V^. 

Divergences in the two-loop contribution V^^ come from G{x, x) and its square. Finite- 
ness of V/3(0, M{(j))) can be shown by first combining and V^^ using the unrenormalized 
form of the gap equation. When the combined expression is written in terms of renormal- 
ized parameters, the remaining divergent integral is cancelled by that of in (3.18). 
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This is another indication that the two-loop contribution must be included for a finite 
self-consistent approximation. The resulting finite expression for V^(0, M(0)) is 

Vp{<j), M{ct>)) = (y° + V'') + V' , (3.19a) 
V' + ^''=2X^- ^M'G{M) - ^4>' , (3.196) 

[A constant term m^/(2A) has been adjusted to obtain (3.19) from (3.12).] However, in 
order to compare the high temperature effective potential with and without the two-loop 
contribution in our later discussion, we still have to extract V^^ from (3.19b). Observing 
that F° must be a function of only, and that in our approximation V^^ does not depend 
on (f) explicitly [since the graph of 0(A) in Fig. 2 does not involve any vertices which depend 
on 0] we obtain, by using the renormalized gap equation. 



Ait y 12 



—G{M)G{M) . (3.20) 
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Clearly the last term in (3.20) is the two-loop contribution. The quantity in the brackets 
is the classical contribution after renormalization is carried out. It is cut-off dependent 
because of the term —X(p^/12, which did not get renormalized due to the structure of the 
gap equation. But the renormalization prescription (3.13) tells us that if Xji is held fixed 
as A — > 00, A approaches 0_. [A necessary condition in a renormalized A0^ theory is 
A < 0.] As shown in the large N studies[17], such theory is intrinsically unstable. On the 
other hand, holding A > implies Ai? ^ as A — > 00. For A > 0, a sensible theory can be 
obtained for a fixed small Aij > as an effective low energy theory, if A is kept fixed at a 
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large but finite value. Such theory requires 

in order to have A > 0, and all momenta, temperature and any other physical mass scale 
must be much smaller than A. 

As shown in Fig. 3, the zero-temperature phase structure of the effective theory with 
finite A is similar to that of perturbation theory: there exists a minimum at a non-zero 
value of 0[16]. We shall consider such an effective theory and take our temperature to be 
T « A. 

C. High Temperature Effective Potential 

Our main interest is in the form of the high temperature effective potential. We first 
find the high temperature gap equation by expanding the integral expression of G{M), i.e. 
the second term in (3.15) at high temperature. Since the basic mass scale in the problem 
is M, we consider an expansion in M^/T^ << 1 



G{M) = 



1 1 M „/M 



r2 



(3.22) 



[We have chosen our renormalization scale /j, to be mR.] Then the high temperature gap 
equation takes the form 



= rfi^ {4>) + -^MT + o( XrM^ In T ) . (3.23) 
24 oTT y J 

Prom the solution of this equation one finds that for a small coupling Xr « 1, the 

condition M'^/T'^ << 1 is consistent with 

m2(0)/T2 « 1 , (3.24) 
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which is exactly the required condition for high temperature expansion of the perturbative 
calculation in Ref . [2] . 

Now we return to Eqn.(3.19) and (3.20). The high temperature expansion of the 
effective potential can be obtained using the high temperature expansion of G{M) and 
also the high temperature expansion of the perturbative one-loop effective potential of 
Ref. [2] by replacing the tree- level effective mass m{(f)) by M 

Vf3{(l>, M{^)) = + + V^^ , (3.25a) 



- 2- 



m 



2 n2 

R 



V 



II 



90 



' 8 



24 



144 247r 



+ 



12^ 

127r 
167r2 



+ 0(M4 In ^), 
+ 0(M^ In — 



(3.256) 

(3.25c) 
(3.25d) 



IV. Discussions and Conclusions 

In order to understand the structure of the effective potential in our approximation, 
we shall first consider the non- linear aspects of the gap equation (3.6), which in the high 
temperature limit can be expressed as 



M\4>) = ^\4>) + - ^TM{cj>) 



(4.1) 



Consequently, M((/)) can be expanded for small \r as 



M(0) = Mi,(0)|l- 



167rMi((/)) 



+ 



\rT Y 



(4.2a) 



where 
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Mz,(0)^ym2(0) + ^T2 (4.26) 

solves the linearized high temperature gap equation, i.e. (4.1) without the last term. 

When the one-loop contribution in (3.25c) is rewritten using the gap equation, we 
obtain 

-4-'^'^- J^M(.).3 - ^ . . (4.3) 

We observe that the term linear in M, namely the third term on the right-hand side 
of Eqn.(4.3), arises from the non-linearity of the gap equation, i.e. from the last term in 
(4.1). If we were to use the linearized gap equation, without this term, the first non-trivial 
correction to the perturbative one-loop effective potential would be given by the term cubic 
in M ((/>). However, at high temperatures the leading non- linear correction is of the same 
order as the term cubic in M(0) in (4.3); in fact from (4.1) we have 

'"■Mm^)/(^)^'^.-Om (4.4) 



19271 //V 1271 J 192 

for T » 0, and such a term could alter the nature of the phase transition. When 
we include the two-loop contribution given in (3.25d) the MT^ term disappears and the 
high temperature effective potential in our approximation is (neglecting (/)-independent 
contributions) 

F(0) = F°(0) + (|^ME(0) - ^M|(0)) (^1 + 0{Xn)^ + 0(^ InT) . (4.5) 

In Fig.4, the effective potential of Eqn.(4.5) is shown as a function of </> for five different 
temperatures close to the critical temperature. It is evident that there is a temperature 
such that there are two degenerate minima. 
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The above analysis of our consistent approximation shows that improving the per- 
turbative one-loop effective potential using the non-linear gap equation clearly leads to 
an erroneous result [18]. Therefore, one must use a self-consistent method which relates 
the effective potential and the gap equation. On the other hand, we also see that, due 
to the cancellation of the leading non-linear effect, one can in fact obtain a consistently 
improved effective potential by improving the perturbative one-loop effective potential us- 
ing an effective mass Ml(0), which is the solution of the linearized gap equation. Such 
improved perturbation theory, using the effective mass squared shifted by a 0-independent 
amount proportional to T^, was first suggested by Weinberg[3] and later further studied 
by others[19]. If the cancellation of the leading non-linear effects in our approximation 
is a general feature, occurring even in gauge theories, the one-loop improved effective 
potential in the standard model first calculated by Carrington[7] would be a consistent 
approximation. We plan to clarify this in a future publication. 
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Figure Captions 
Fig. 1: Examples of daisy (a) and superdaisy (b) graphs. 

Fig. 2: The two- and three- loop graphs contributing to F^ {(f), G). They are two-particle 
irreducible and their lines represent the full propagator G. 

Fig. 3: Vji{(p,M[(p))/mji as a function of 4>/mR at T = 0. Vp illustrated in figure cor- 
responds to \r = 0.05 and ln(A^/m|j) = IGtt^. becomes imaginary for small 
(t)/mR. 

Fig. 4: The high temperature effective potential in our approximation (Eqn.(4.5)) for five 
different temperatures close to the critical temperature. [Here we have chosen 
Xr = 0.05 and ln{A'^ /m%) = IGtt^.] In order to compare the shapes between 
different temperatures, we shifted the T/g's by 0-independent amounts so that 
= for (p/niR — 0. From the figure it is evident that at a certain temperature 
there are two degenerate minima. 
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